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SECTION 1 


INTRODUCTION 

In converting solar radiation to electric power in space, there is 
a high premium on the ability of the converter to operate at high tem- 
perature. There are two fundamental reasons for this. First, higher 
operating temperatures permit using a higher concentration ratio for 
the incident solar radiation. This favorably affects the cost and 
weight of large systems since the cost and weight of the concentrator 
per unit surface area are generally less than those of the collector. 
Second, operating the converter at high temperature permits using 
smaller radiators for rejecting the waste heat. (The radiator is 
typically the heaviest component of projected space power systems.) 

These premiums on high temperature apply to solar electric converters 
in general. In the heat-engine approach to solar electric conversion, 
there is the additional advantage that the working cycle efficiency is 
higher at higher temperatures. 

This report presents the results of a two-year study aimed at 
developing an ultra-high-temperature solar electric heat engine converter. 
In the converter concept, an alkali plasma serves as both an efficient 
high-temperature collector of solar radiation as well as the working fluid 
for a high- temperature working cycle. The working cycle is a simple 
magneto-hydrodynamic (MHD) Rankine cycle employing a solid electrode 
Faraday MUD channel. 

The simplest ana most attractive solar-powered MHD cycle for high- 

power (£10 kW) solar electric converters in space is a plasma-phase 

MHD converter utilizing a pure alkali Rankine cycle. Solar radiation 

concentrated by a factor of several thousand would heat cesium or 

3 

potassium to Cj 3000°K at a pressure of about 10 Torr. The hot plasma 
would flow through a nozzle and across an VLO-kG magnetic field, thereby 
driving an electrical current transverse to the flow direction. At the 
above temperature and pressure, a pure cesium plasma will have a high 
enough electrical conductivity to enable electric power to be extracted 
efficiently from the plasma with a simple Faraday MHD channel employing 
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solid electrodes."* - The plasma would remain sufficiently conductive to 

extract power in this way down to a temperature of 'VL800°K. The net 

cycle efficiency is expected to be ^15%, and, at a 1000°K condensor 

temperature, the reject heat can be radiated away at a rate per unit 

2 

radiator area of %60 kW/m . 

An alternative to radiating away the reject heat from the plasma 

MHD cycle would be to use the reject heat to drive a liquid-metal MHD 

bottoming cycle. The projected efficiency of one liquid-metal MHD cycle 

2 

operating between 300 and 1000°K is calculated to be ^60%. Such a 

binary cycle could retain the simplicity, ruggedness, and relative 

inexpensiveness of MHD and would operate with an overall efficiency of 

^66%, but it would do so at the expense of requiring a considerably 

higher radiator area in spite of the higher efficiency because of the 

fourth-power dependence of the radiation rate on temperature. 

The above approach to space solar power was originated independently 

3 

by Hughes Research Laboratories (HRL) and by the University of 
Washington's Aerospace and Energetics Center about two years ago. We 
have achieved during those two years three major milestones towards the 
realization of this type of converter. First, we formulated a complete 
and versatile theoretical model of the physics of coupling solar radia- 
tion to a cesium plasma. The model, in the form of a computer program, 
treats a variety of experimental conditions and is easily adapted to 
treat the other alkalies. Second, we have considered in some detail 
three MHD working cycles.; plasma MHD using a single alkali species 

cycle; two-phase liquid-metal MHD using either a single or two alkali 
2 

species cycle; and an alkali rare gas (Cs-Xe) plasma, MHD cycle in 
which excimer laser power is extracted in place of electric power. ^ 
Finally, we recently demonstrated heating of a cesium plasma in a labora- 
tory experiment using a solar simulator light source. ^ The measured 
temperature of the plasma was within 20% of that predicted by our 
theoretical model for the conditions of the experiment. 

The major theoretical and experimental research milestones attained 
to date on this program are presented in Sections 2 and 3. Important 
research objectives remain for establishing firmly the feasibility of 
this concept for space solar power, and these are discussed in Section 4. 
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SECTION 2 


THEORETICAL PROGRAM 


A. IDENTIFICATION OF HIGH-TEMPERATURE COLLECTOR 

To achieve highest efficiency in the heat-engine approach to 
solar-electric conversion, it is advantageous to deposit the solar 
radiation into the working fluid at as high a temperature as possible. 
For this temperature to approach the upper limit of the sun’s tempera- 
ture, the radiation must be absorbed directly in the working fluid in a 
plasma state. Therefore, a plasma in the temperature range of 3000 to 
6000°K needs to be identified that presents a broadband absorption cross 
section throughout the solar spectrum that is large enough to permit 
using moderate densities and dimensional scales for efficient coupling 
of the radiation to the plasma. 

The primary sources of broadband absorption of visible radiation in 
plasmas are inverse bremsstrahlung (free-free electron transitions) and 
photo- ionization (bound-free electron transitions). If one assumes 
thermal ionization, as described by Saha's equation, then the absorption 
cross section (per ground state) caused by free-free transitions in a 

g 

hydrogen-like plasma is 


a 


ff 




exp (-V x /T) 
2 (V ] ./T) 


2 

cm 


( 1 ) 


where is the ionization potential, T is the plasma temperature, and 
hv is the photon energy. At a temperature near 5000°K, even with a low 
ionization potential species, the cross section is quite small; there- 
fore, very high densities and associated high plasma pressures will be 
needed to obtain a practical value for the free-free absorption coeffi- 
cient in the solar spectrum (hv ^ 1 to 3 eV) . 

On the other hand, the absorption cross section caused by bound- 
free transitions can be much larger in the visible than the free-free 
transition cross section. Combining Kramer's formula for the hydrogenic 
bound-free cross section per excited state with an assumed Boltzmann 
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distribution for the excited state population yields the cross section 

, 8 
per ground state as 



where n*(i) is the effective principal quantum number, L(i) is the total 
angular momentum quantum number, and E(i) is the binding energy of the 
i C ^ excited state. 

For the temperature of interest (3000 to 5000°K) , cesium exhibits 
the largest average bound-free cross section through the solar spectrum. 
Cesium also has other important continuum absorption bands iu the red 
portion of the solar spectrum as a result of cesium dimers. These 

O 

absorption, bands extend from about 6000 A to 1 ym, which, conveniently, 
is where the bound-free cross section is falling off. The aimer absorp- 
tion cross section may be computed from the quasi-static theory of line 
a 

broadening. This technique, which requires knowing the potential energy 
curves for the transition, partitions the molecular states into a set of 
substates characterized by their internuclear separation and populated 
relative to their dissociation products according to a classical conical 
ensemble. Thus, the concentration of dimers in the ground electronic 
state with an internuclear separation between R and R + dR is"^ 

d[X] = 4 R 2 dr(g x /g f ) exp(-WX(R)/kT)[Cs]*- , (3) 

where [ ] indicate species concentration; g and g f are the degeneracies 
of the molecular X state and its parent atomic state, respectively; 

WX(R) is the energy of the molecular substate relative to that of its 
parent atomic state; k is Boltzmann's constant; and T is the gas 
temperature. 

The absorption coefficient caused by a ground-state dimer transition 
9,10 

can then be wrxtten as 

B Di = 2 iik (-WX(R) /kT) cg x /g f ) [CS ] 2 , (4) 
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where A is the Einstein coefficient for the transition (assumed 
independent of R) , and \ and V are the wavelength and frequency of the 
radiation. 

Rather than computing the dimer absorption by simply inputting the 
measured values of the dimer absorption cross section at an arbitrary 
temperature (as done, for example, in Ref. 11), we chose to compute the 
temperature-dependent absorption cross sections according to a four- 
Morse-potential state model of the cesium dimer molecule. Not only does 
this allow accounting for the temperature dependence of the cross section, 
it also assures that re-radiation is accounted for properly on the dimer 
transitions. The 8,000 to 10,000 A X-A absorption band can be modeled 
fairly exactly in this way since it is due to transition to a single 
state. And although the shorter wavelength bands are actually caused by 
several overlapping transitions, the absorption can be modeled to an 
accuracy sufficient for our purposes by constructing three hypothetical 
Morse potential states with the parameters given in Table 1. Our com- 
puted absorption cross sections (shown in Figures 1 and 2) agree well 

13 

with Weschler's measured values. Line absorption is not included in 
the theoretical model since it contributes negligibly to the net energy 
flow in the plasma. 


Table 1. Morse Potential Parameters for the Four-State 
Model of CS2 Absorption Spectra (X-state 
parameters from Ref. 12) 


State 

_ -1 
D, cm 

R , A 
e’ 

W , cm ^ 
e’ 

* 

Atomic State 

X 

3600 

4.47 

42 

6S 

A . 

5100 

5.25 

34 

6P 1/2 

B 

1510 

4.47 

34.2 

6P 1/2 

G 

500 

5.7 

15 

6P 3/2 
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4 STATE MORSE POTENTIAL i 

MODEL I 

WESCHLER'S MEASUREMENTS I 

I 




0.40 0.44 0.48 0.52 0.56 0.60 0.64 0.68 0.72 0.: 

WAVELENGTH, ym 

Figure 1. Computed and measured absorption cross sections of cesium 
vapor: it- 40 to 0.75 ym (measured values are from Ref. 13). 




4 STATE MORSE POTENTIAL 
MODEL 

WESCHLER'S MEASUREMENTS 


0.74 0.76 0.78 0.80 0.82 0.84 0.86 0.88 0.90 0.92 0.94 0.96 0.98 1.00 

WAVELENGTH, ym 

Figure 2. Computed and measured absorption cross sections of cesium 
vapor: 0.74 to 1.00 ym (measured values are from Ref. 13). 
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The net absorption coefficient due to both the bound- free transitions 
and the dime? transition illustrated in Figure 3 in an LTE cesium plasma 
at a temperature of 4000°K and a pressure of ^-2 atm is shown in Figure 4, 
where it is compared with the solar spectrum. Both the spectral profile 
and the magnitude of the cross section are shown to be well suited for 
the use of cesium at moderate densities as an efficient sol's? radiation 
receiver. These features and its relatively high vapor pressure at prac- 
tical operating temperatures, its high electrical conductivity, and its 
availability identify cesium as an excellent working fluid for an 
advanced, high- temperature solar-electric heat engine. 



CESIUM ABSORPTION SOLAR FLUX. mW/(cm 2 • 100 A) 

COEFFICIENT, cm -1 


6810-4R4 




Figure 4. Solar spectrum and cesium absorption spectra. 
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3. 


IDENTIFICATION OF A HIGH-TEMPERATURE WORKING CYCLE 


Having identified a cesium plasma as the optimal high- temperature 
collector for an advanced solar heat engine, our next task was to iden- 
tify an optimal working cycle for converting the plasma thermal energy 
to electric power. The simplest and most direct way is to allow the 
plasma electrons themselves to carry the load current. There are only 
two mechanisms by which charge separation leading to net current flow 
out of the plasma can be accomplished in a thermal plasma: 

• Charge separation caused by ambipolar diffusion 
(Weymouth plasma diode, TELEC-llke devices) 14, 15 

® Charge separation caused by JxB forces (MHD) . 

* 

In the first method, electric power coupled out of the plasma via 
ambipolar diffusion is essentially a surface rate that is independent of 
density and depends inversely on the smallest dimensional scale of the 
plasma: 


P 


amb 



fD (1.4 • 10" 23 ) T 
a e 

x 


(5) 


where f is the degree of ionization of the plasma, D is the ambipolar 

-1 a 

diffusion coefficient per unit density in cm /sec, T^ is the electron 
temperature in °K, and x is the small dimension of the plasma in cm. ' 

We estimated that, at the temperatures of interest (^3500°K) , the 
cesium plasma would be losing energy to radiation at an effective sur- 
face rate of about 10% of the blackbody rate: 


P 


rad 



0.1 (5.7 


10“ 12 ) T 4 


( 6 ) 


At T m 3500°K, D is on 
a 

for ambipolar diffusion 
electrical power out of 


the order of 10 2 ^ cm Vsec. 1 ^ Thus, at 3500°K, 
to compete with radiation in efficiently coupling 
the plasma, we must have 


15 


(7) 


J/x £ 10 2 cm -1 

Unfortunately, this condition is totally incompatible with efficient 
coupling of the sun's radiation into the plasma. Therefore, ambipolar 
diffusion cannot be anticipated to be an acceptable thermal-to-electric 
conversion mechanism when applied to a solar-sustained plasma. This 
leaves MHD as the apparent optimal working cycle. 

In the MHD cycle, the plasma thermal energy is converted to trans- 
lational energy by expansion of the plasma into a lower pressure region 
as in a turbine cycle but with the conversion of the plasma flow energy 

- 4 - - 4 * 

to electricity taking place by action of the JxB force on the plasma 
rather than by a turbine. In this way, power can be coupled out of the 
plasma at a rate that can be made to compete with the radiation losses 
from the plasma. 

The most important plasma parameter for an MHD working fluid is its 

conductivity. The electrical conductivity of a 1 atm pressure cesium 

plasma in the temperature range from ^1800 to 3000°K is on the order of 
2 3 17 

10 to 10 mho/m (see Figure 5). At a magnetic field of ^10 kG, plasma 
conductivities of this order will require MHD duct lengths of ad m for 
significant extraction of electrical power from the plasma enthalpy (see 
Figure 6) 

The Hall parameter is the other important parameter for an MHD 
plasma (electron-cyclotron frequency times mean-free time) . The Hall 
parameter for a 1 atm pressure solar-sustained cesium plasma at tempera- 
tures of a,3000°K is much greater than 1 for an ad0-kG field. This per- 
mits using a Faraday generator configuration, which does not require the 

use of segmented electrodes. Enthalpy extraction efficiencies of 15% 

1 

have recently been demonstrated with this type of generator. 

The optimal thermodynamic MHD working cycle was identified to be a 
simple Rankine super-heat cycle. A conceptual description of the cycle 
is presented below. 

Solar radiation is admitted through a sapphire window and brought 
to a focus in a region adjacent to a reservoir of metallic cesium. The 
proximity of the cesium reservoir to the focal region is chosen such 
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dc ELECTRICAL CONDUCTIVITY, mho/m 




! I vju i *.■ 6. Required MHl) duct length for conversion of 
pliisnui cn tlnil pv to electric power from 
Kef, Id. 
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Chat: radiation from the plasma will provide sufficient heating of the 
reservoir to give the desired vapor pressure. The sapphire window is 
placed far enough away from the focal region and cesium reservoir so 
that the condensation on the window can be prevented by local window 
heating. 

Within the focal volume, the cesium vapor is brought to a tempera- 
ture of about 3000°K by absorption on dimer transitions and excited 
state photoionizing transitions. The focal region forms the entrance to 
a plasma MHD duct. The duct is profiled so that the plasma will expand 
through it at a roughly constant velocity. The magnetic field strength 
within the duct is on the order of 10 kG which is provided by a super- 
conducting magnet. It should be possible to control cesium condensation 
on the insulating walls of the duct by using either local heating or a 
buffer gas boundary layer. 

After passing through the MHD duct, the cesium vapor is recondensed 
at a temperature of about 1000°K, corresponding to a pressure difference 
through the generator on the order of 1 atm. The molten cesium is then 
returned to the focal region reservoir by an electromagnetic liquid- 
metal pump. 

A rough estimate of the scalability of the concept was obtained 
from the constraint that the aspect ratio of the MHD duct dimensions not 
exceed a value that would lead to excessive conduction losses. A pre- 
viously used criterion would require that the ratio of the cross-sectional 

18 

dimension to the length of the MHD duct not be less than M./20. The 
minimum MHD duct length required for extraction of all of the plasma 
enthalpy is set by the stagnation pressure, the plasma conductivity, 
and the magnetic field. The stagnation pressure, of course, also 
determines the deposition length of the radiation in the plasma. Assuming 
a solar radiation concentration ratio of 10 (yielding a plasma tempera- 
ture of ^3000°K, which corresponds at 1 atm to a plasma conductivity of 
3 17 

'VlO mho/m) and assuming a 10-lcG field, the minimum duct length 

18 

appears to be *V1 m, corresponding to a minimum aperture of 't'lQ cm or a 
minimum input power level of n>io kW. 
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C. RADIATIVE ENERGY BALANCE RESULTS 

A. theoretical radiation energy balance model was developed to 
investigate the limitations imposed by radiation losses alone in achieving 
high temperatures in a solar-sustained cesium plasma. The model was 
later expanded to include thermal conduction, convection, and gas flow, 
as would occur in the working fluid of a solar electric heat engine 
(see Section 2.D) . 

The basic assumptions of the model are that local thermodynamic 

equilibrium exists and that the electron and heavy particle temperatures 

19 —3 

are equal. Since, at the cesium vapor densities of interest (^10 cm ), 
electron/heavy-particle collisional energy transfer rates will dominate 
radiative and other convective energy transfer rates, these assumptions 
will be vaoid; also, reabsorption of radiation is neglected. 

In the static plasma model, the plasma temperature is assumed to be 
independent of position in the plasma. The temperature is computed by 
balancing the energy absorbed on the three dimer transitions and eight 
photo-ionization transitions of cesium against energy lost to spontane- 
ous emission and radiative recombination into these same states. The 
radiative rates are calculated on the basis of detailed balancing applied 
to the absorption cross sections. Rates into and out of states higher 
than the eighth excited state of cesium were verified to contribute 
negligibly to the overall energy balance in the plasma. 

The energy balance calculations were carried out with a computer 

program. The results are shown in Figure 7 as a plot of the required 

o 

solar flux (in W/cm“) for maintaining a cesium plasma against radiative 

losses versus the attainable plasma temperature. The solar spectrum was 

taken as the spectrum of a blackbody at a temperature of 5800°K. The 

2 

concentration ratio required to achieve these fluxes for a 1 kW/m solar 
flux assumed incident at the earth is shown on the right ordinate. The 
maintenance flux for a perfect blackbody is shown for comparison. For 
these results to be valid for the static plasma, the energy loss to 
thermal conduction must be small compared with the radiative losses, and 
the spatial extent of the plasma along the radiation propagation direc- 
tion must be small compared with the smallest penetration depth of the 
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solar spectrum in the plasma. These conditions will be approximately 
satisfied for a 1-cm-diameter spherical cesium plasma at the assumed 
density for the temperature range from ^2500 to 4000°K. As shown in 
Figure 7, higher plasma temperatures result at lower oven temperatures 
in the purely radiative energy balance case because the lower fractional 
concentration of dimers causes the spectral range over which most of the 
absorption occurs to move to the shorter wavelengths associated with the 
bound-free transitions. 

It is possible to understand the general form of the result pre- 
sented in Figure 7 on the basis of a simple physical argument based on 
the principle of detailed balancing. The argument goes as follows. 

Consider a system that is absorbing and emitting radiation on one 
transition only. If the system is in steady state, then the energy 
balance condition is 


P(v)a(n^ - n u ) = n u A(hv) 


( 8 ) 


where n u and n^ are the population densities of the upper and lower level 
for the transition, O is the absorption cross section for the transition, 
A is the spontaneous radiative rate, and P(v) is the input radiative 
power flux to the system at the transition wavelength. If the level 
populations are in thermal equilibrium at a temperature T, then, from 
detailed balancing, P (v) is the Planck spectrum: 


P T (v) = 


3 

2hv 

2 

c 


(e hWkT . ir l 


(9) 


Alternatively, if blackbody radiation from a radiator of temperature T 

K 

is to sustain the system at a temperature T, then the energy balance 
condition is 


CP„ (V) 
R 



n 

u 



= P T (v) 


( 10 ) 
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where C is a factor which is equal to one at the surface of the radiator. 
Obviously C = 1 for T = T R . If T R / T, then 

P T (v) 

c = JTJy) ‘ (11) 

R 

From the form of the blackbody spectrum (Eq. 9) , C < 1 for T < T R and 
C > 1 (optically impossible) for T > T R . The above argument continues 
to hold (term by term) for a system absorbing blackbody radiation on 
more than one transition and, in particular, on the continuum bound-free 
and dimer transitions. The essential condition for the argument to be 
valid is that the system energy balance be dominated by radiation, that 
the input radiation be in a blackbody spectrum, and that the system be 
in thermal equilibrium at some temperature T < T R . 

In the case of several transitions or continuum transitions, the 
value for C will be determined from a spectrum integration of Eq. 11. In 
this case, one can write approximately: 

( exp (hv/T D ) - l) 

C^< \ , (12) 

(exp(hv/T) - 1) 

where V is an average frequency of the absorption spectrum. This is 
the general form of the solid curves in Figure 7. 

The re-radiation spectrum by the continuum absorption transitions is 
shown in Figure 8 for a selected plasma temperature. Although this 
radiation dominates line radiation in controlling the plasma temperature, 
line radiation will be a more useful signature of the plasma temperature 
due to its higher spectral brightness (see Section 3.C). 
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RECOMBINATION RADIATION 
STRENGTH (REL. UNITS) 


8547-2 



Figure 8. Theoretical continuum re-radiation spectrum 
(minus line emission) . 
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DIMER RADIATION STRENGTH 
{REL. UNITS) 




D. TOTAL ENERGY BALANCE RESULTS 


Lower oven temperatures result in reduced radiative energy deposition 

in the plasma, which will eventually make losses other than radiation 

dominate the energy loss from the plasma. Therefore, we have added both 

thermal conduction and free convection losses to the model. For the 

densities, convection velocities, and dimensional scales of the present 

experiments, the ‘estimated Reynolds number for the plasma is quite small 

(•A>100) , permitting the use of a kinematic value for the thermal con- 

ductivity coefficient for cesium vapor (VL5.1 • 10 W/cra°K). Since 

the free convection flow was not viscosity limited for the dimensional 

scale of the present experiments, it was computed from buoyancy forces 
20 

alone. These nonradiative losses control the dependence of the plasma 
temperature on vapor pressure in our experiments for oven temperatures 
below about 800°K, as shown in Figure 9. 

The model also accounts for the convergence of the incident radia- 
tion and computes the plasma temperature as a function of the propagation 
distance of the radiation into the vapor. The results, shown in Figure 10, 
show, for example, that higher plasma temperatures are attainable by 
positioning the focal volume as close as possible to the front boundary 
of the vapor. 

Finally, for treating the flow of the plasma into a working cycle, 
we added a one-dimensional flow of the plasma along the propagation axis 
of the incident radiation to the model. The results are shown in 
Figure 11, in which it is assumed that the solar concentration ratio is 
3000. Figure 11 shows the temperature, the fraction of incident radia- 
tion power remaining in the beam (transmitted power) , and the fraction 
of incident radiation power contained in the plasma enthalpy flow 
(coupling efficiency) . These are shown as functions of z for chosen 
values for the input flow velocity and ambient cesium vapor pressure. 

The coupling efficiency equals one minus the sum of the fractional 
power transmitted and the fractional power radiated. The results show 
coupling efficiency to be maximum at the same downstream location as is 
temperature. In a heat-engine application, this location would be the 
optimal point at which to extract the cesium plasma working fluid from 
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8547-3 



Figure 10. Plasma temperature as a function of distance from the front 
boundary of the vapor for two focal spot positions (theory). 
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Figure 11. Plasma temperature, coupling efficiency, and transmission 
coefficient as a function of propagation distance in a 
cesium plasma from the one-dimensional convective model. 
(Oven temperature = 1000°C a solar concentration ratio = 
3000.) 
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the radiation interaction zone. This interaction length scales 
proportionally with the flow speed and somewhat faster than inversely 
with the density (because of the influence of the dimer). There is a 
well-defined tradeoff between peak coupling efficiency and peak plasma 
temperature: higher flow velocities yield higher coupling efficiencies 

but (as expected) lower peak temperatures. Since, in a heat-engine 
application, the working cycle efficiency will generally increase with 
the input temperature of the working fluid, the requirement for maximum 
overall efficiency will determine an optimal flow speed through the 
radiation interaction zone. 

All of the above plasma maintenance calculations are done with a 
computer program written for a PDP-10 computer. The flow chart 
for the computer program is presented in Appendix A. 

E. SOLAR PUMPED LASER RESULTS 

The concept of a solar-pumped MD laser was first disclosed at the 

1978 NASA Ames Conference on Radiation Energy Conversion and in the 

21 

HRL technical proposal for this contract. It was pointed out then that 
a population inversion on the cesium dimer A-X laser transition will 
develop in a solar-sustained cesium plasma if the electron temperature 
can be made to exceed the gas temperature by at least the fractional red 
shift of the transition. It was suggested that such a decoupling of the 
electron and gas temperatures could be made to occur in a nozzle flow of 
the plasma across a magnetic field. Thus, a closed Rankine cycle solar- 
pumped laser was identified that is similar in concept to the MHD solar- 
electric converter with laser power replacing or supplementing electric 
power as the work output of the cycle. The laser characteristics for 
this type of transition would include^ 

o 

• High saturation power (^50 kW/cm ) 

• Large, high-power system capability (low gain, fewer 
parasitic problems) 

• Broadband tunability (AA 'v 1000 A) . 
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During this program, we quantified these speculations with a 
mathematical model and expanded the lasing transitions considered to 
include the CsXe A-X excimer band. The excimer laser transition will 
offer essentially the same advantages as listed above for the dimer 
transition. The results of the model appear favorable for the realiza- 
tion of laser action on the excimer transition in this system. 

An energy level diagram showing the CsXe excimer and Cs 9 dimer laser 
transition is shown in Figure 12. The fractional red shift is ^0.07 
for transitions from the bottom of the excimer A state and M).23 for 
transitions from the bottom of the dimer A state. The two- temperature 
equilibrium model would predict a population inversion to occur on the 

excimer band for T /T £ 14 and on the dimer band for T /T £ 4.3. 

eg e' g 

These temperature ratios are easily attainable in a supersonic MHD flow 

(as will be seen below). However, for the plasma conditions we consider, 

several factors will conspire to invalidate the use of the simple two- 

temperature equilibrium model to determine the net gain on these laser 

transitions. These factors are discussed below. 

The starting point for calculating the electron and gas temperature 

in a plasma flow through a short-circuited Faraday MHD channel is an 

23 24 

electron energy balance equation derived by B. Zanderer et al. * and 
utilized in their program to generate 10-um laser action in an MHD 
channel: 


T 

' Y- - 1 = f M 2 y^d" 1 , (13) 

g 

where y is the ratio of specific heats; 3 is the Hall parameter; M is the 

Mach number of the flow; and 6 is an inelastic collision energy loss 

factor, which can be taken as equal to one plus the ratio of the inelastic 

electron collisional energy loss rate to the elastic collisional energy 
23 24 

loss rate. ’ If a Maxwell-Boltzmann energy distribution is assumed 

for the electrons and the excited state populations of cesium, then 6 

25 

may be estimated as 
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Figure 12. CsXe excimer and Cs 0 dimer laser transitions. 
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(14) 


■V “ 1 + 


(6.38 • 10 ? )(2 T e ) 3/2 expC-E^T^Cl + E^T^) 

V 2 W Ui '(W®. - v dE 


where is the energy of the first excited state of xenon; is its 
excitation cross section; is the electron, elastic collision cross 
section; c is the speed of light; and the energies, temperatures, and 
masses are all expressed in eV. 

If we add to Eqs. 13 and 14 the isentropic expansion equation for 

. „ „ 24 

the gas temperature, 

T g /T o - jl + ~ (Y “ 1)M 2 | , U5) 


where T is the stagnation temperature, then we may solve for the electron 
0 

temperature and gas temperature as a function of the Mach number and the 


We use Y " 5/3, 0 


=17 2 “1 

10 cm eV \ 


el 


10 ^ cm 2 , and da^/dE 


Ei 


The collision cross sections <J el and dc^/dEl^ 


Hall parameter. 

* 3 

were estimated from data on lithium in Ref. 26. The results for the 

electron temperature and the gas temperature are shown in Figure 13 as a 

function of Mach number for 3 “ 10 and T q = 0.2 eV. 

The small-signal gain coefficients on the excimer and dimer bands 

10 

can be computed from the quasi-static theory of line broadening: 


gain 


,2 - r - 1 

A T 
8tt 


I dfA] _ £a d [X] 1 dR 
dR ” g,. dR l dv 


(16) 


’X 


where A and T are, respectively, the wavelength and radiative lifetime 
of the transition; V is the frequency of the transition at the inter- 
nuclear separation R; [A] and [X] symbolize the concentration of the A 
and X states, respectively; and g A and g x are their degeneracies. The 
theory assumes that the molecular states are populated as a function of 
R according to a cononical ensemble so that we can write 
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Figure 13. Electron and gas temperature as a function of Mach 
number for a cesium-xenon plasma flow in Faraday 
MHD channel. 



32 



gai.nO) “ i(\) 


f-U) «A [X! «P(-V x (X>/«.> ) 

\h (g A,f /g A^ s x k x (g X,f /g X J ) 


( 1 ?) 


where and are Che A stace and X scare equilibrium constants, the 
gj^s are the degeneracies of the atomic parent states, V^(A) and V^(A) 
are the potential energies of the A and X states at the internuclear 
separation corresponding to the transition wavelength A, k is Boltzmann's 
constant, and >t(\) is the stimulated emission coefficient given by 

“<« • (sfe) r2 f • <“> 

At the temperatures and densities of interest here, the population 
of the CsXe A state is controlled by the heavy particle dissociation 
rote. In this case, both the A state population and the X state popu- 
lation will be in equilibrium with their parent atomic state, making the 
two- temperature equilibrium model for the population of the upper and 
lower molecular states of the excimer transition valid. However, in 
determining the net gain on the excimer level, one must also include the 
absorption on the X-*-A dimer transition since this absorption band over- 
laps the excimer gain band. Under these conditions, we write the 
excimer gain as 


gain ex CX) 


a ex (X) 


j[Cs*][Xe] 


exp(-V A (A)/kT R ) 

(g A,f /g A^ 


S A 

“TT [Cs][XeJ 
8 X 


exp(-V (A)/kT ) 
- (s x >s x ) 





[Cs 2 ] 


3/v «p(-V x /kT ) 


K X 8 X ^X.f'V 


'di 


( 19 ) 


where [Cs*] represents the concentration of cesium atoms in the 6P reson- 
ance state. 
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On the other hand, the Cs 2 A state population will not in general 
b? in thermal equilibrium with the parent atomic state at the tempera- 
tures and densities of interest. The relatively large dissociation 
energy of this state will cause the collisional dissociation rate of the 
state to fall below its radiative rate at the lower temperatures and 
densities associated with high Mach numbers. With the radiative rate 
contributing to the net loss rate of the steady-state dimer A state popu- 
lation, Eq. 17 for the gain coefficient becomes 


gain di (A) = a di (A) 


(r A [Cs*J[Cs][Xe] 

exp(-V A (X)/kT g ) 

( 8 A,f /s A 

K.CT 1 + F. [Xe]/K A ) 


_ g A fCs] 


3 X s X,f' e X 


ex P (-V x (A)/kT g ) 


'di 


( 20 ) 


where F. is the three-body association rate coefficient for the formation 
of the A state. The scaling arguments presented in Ref. 27 lead to an 
estimate of F^ - lx 10~ 2 ' 1 ‘cm f ’sec~' 1 ‘ and 8 x 10 -22 exp (-V A /kT g ) , where 

is the potential energy at the bottom of the dimer A state well. 

For the conditions we consider, the dissociation time of the dimer 
ground state will be comparable to the flow time of the plasma through 
centimeter-scale zones. Thus, the population of the dimer X state will 
not always be in equilibrium with the monomer concentration and in gen- 
eral will be given by 

ICs 2 ] = K x [Gs] 2 {1 - exp (F diss [Xe] z/V z )} , (21) 


where V and z are the flow velocity and distance along the flow, respec- 
z 

tively, and r ddsg is the dissociation rate constant for the X state and 
is related to the recombination rate constant and equilibrium constant 
by F 


= K v F . 
diss X X 


From Refs. 27 and 28, we estimate F 


From Ref. 13 we find K x = 2.7 x 10 


-23 


X 


n in -32 -6 -1 

1 x 10 cm sec 


exp(-V x /kT g ) . 
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The maximum gain coefficient on the excimer and dimer bands occurs 

near 0.96 um and 1.2 ym, respectively. Inserting in Eqs. 18 through 20 

o 

the R“ and dR/dv values at these positions on the potential curves and 
the lifetime and relevant degeneracy factors from Ref. 10, we obtain 

gain ex (0.96 ym) - 1.5 • 10“ 37 {exp(0. 05/T g ) exp (-1. 45/T e ) 

- exp(-0.1/T ) } [CsJ [XeJ 

O 

- 2.0 • 10 {0. 25/(2, 7 • 10 )}[Cs 2 ] (22) 

gain^^(1.2 ym) * 2.0 * 10 


0.25 • 1 ♦ 10" 31 exp(-l. 45/Te) 


(8 * 10 23 ) 5.9 * 10 7 + 1 • 10" 31 [Xe] / (8 • 10 -23 exp (0, 63/Tg) ) ] 


[CsT [Xe] 


• 2.0 


10 -37 [Cs ] 

((2.7 * 10“ )) 4 


(23) 


We have assumed that the population of the 6P resonance level at 1.45 eV 

is in equilibrium with the monomer ground state at the electron tempera- 

—3 

ture. All temperatures are in eV, concentrations are in cm , and the 
gain coefficients are in cm™\ 

The stagnation cesium density for the plasma flow was chosen to be 
compatible with the coupling of solar radiation to the plasma through 
dimensional scales on the order of several centimeters, while the xenon 
density was chosen to maximize the excimer gain for minimum magnetic 
field and Mach number. The chosen stagnation densities were: 

[Cs] D = 3 • 10 17 cnf 3 (24) 

on _o 

[Xe3 o « 1 * 1(T U cm J . (25) 
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At an assumed temperature of 0.2 eV for the solar-heated plasma, these 
densities correspond to a stagnation pressure of *^20 atm. 

Figure 14 shows the results for the small-signal gain coefficients 
on the dimer and excimer bands computed for these conditions as a func- 
tion of Mach number using Eqs. 13 through 23. The figure also shows the 
approximate magnetic field required to maintain the assumed value of 10 
for the Hall parameter. These results show that positive gain coeffi- 
cients are available on the excimer band near a Mach-2 flow and above 

about a Mach- 7 flow. The dimer gain remains negative throughout the 
29 

entire range. The minimum gain region around Mach 4 roughly divides 
the kinetic regime for the ground state dimers into a regime where the 
dimers are in equilibrium with the monomers for smaller Mach numbers and 
a nonequilibrium regime for larger Mach numbers. This is illustrated 
further in Figure 15. The Mach-2 excimer gain is shown to be constant 
with distance in the flow, indicating that equilibrium concentrations of 
all species involved in determining the excimer gain have been reached. 

On the other hand, the Mach-7 excimer gain is seen to fall linearly with 
distance, which reflects the increasing concentration of absorbing ground 
state dimers with distance as they approach their equilibrium values at 
some distance much further down stream. 

The high-Mach-number regime for positive excimer gain extends over a 

broader range of flow parameters and does not require the high magnetic 

fields that the low-Mach-number regime does. Also, it represents a 

kinetics regime for the alkali rare gas excimer laser that has not yet 

30 

been explored experimentally. Thus, if these high Mach numbers can be 
attained in this type of plasma, this regime would appear to offer the 
greatest potential for the first demonstration of an alkali rare gas 
excimer laser. 

Arc or chemical heating might be used as an alternative to solar 
heating the plasma. Shock tube experiments utilizing an arc-heated 
plasma probably represent the most expedient means for performing the 
initial experimental investigations of the concept. 
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Figure 14. Small signal gain coefficients and required magnetic 
field versus Mach number for a cesium-xenon plasma in 
a Faraday MHD channel. 
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SECTION 3 


EXPERIMENTAL PROGRAM 


A, RADIATION SOURCES 

Two alternative radiation sources were used on this program. A 
15-in. Fresnel lens solar concentrator mounted on an equatorial tracking 
mount was made available for outdoors use. The equatorial mount was 
large enough to accommodate both the cesium absorption cell and a 1/4-m 
Jarrel ash monochromator and vidicon tube assembly. This apparatus is 
shorn in Figures 16 and 17. The parameters of the Fresnel lens concen- 
trator are presented in Table 2. 

» 

For indoor use, a 1600-W xenon arc lamp adapted with a 10-ln.- 

diameter elliptical reflector for simulating concentrated sunlight was 

purchased by HR as capital equipment for use on this program. This solar 

simulator is shown in Figure 18. The radiation flux profile provided at 

the focal plane of this source is shown in Figure 19. The peak radiation 
2 

flux of ^30 W/cm is equivalent to a solar concentration ratio of about 
300. 


Table 2. Parameters of Solar Concentrator 


LENS DIAMETER 

38 cm (15 in.) 

FOCAL LENGTH 

61 cm (24 in.) 

f 

1.6 

AREA 

1140 cm 2 

TRANSMISSION 

85% 

SOLAR IMAGE DIAMETER AT FOCAL PLANE 

tan 0.5° x 61 cm = 0.53 cm 

SOLAR IMAGE AREA AT FOCAL PLANE 

0.222 cm 2 

MAXIMUM CONCENTRATION RATIO AT 
FOCAL P LANE 

1140 -0.85 = 4365 
0.222 

MAXIMUM COLLECTED POWER DENSITY AT 
FOCAL PLANE (AIR MASS ONE 86.5 mW/cm 2 ] 

4365 > 86.5 • 10'" 3 = 377 W/cm 2 

POWER COLLECTED FOR 1 cm 2 
ROUND CELL (DEFOCUSSED) 

1140 • 86.5 • 10“ 3 • 0.85 = 84 W 

POWER COLLECTED FOR 1 cm 2 

SQUARE CELL (DEFOCUSSED AND MASKED) 

84 W • = 33 W 

POWER DENSITY CONCENTRATION RATIOS 

1 cm 2 ROUND CELL 1140 • 0.85 = 969 
1cm 2 SQUARE CELL 969 = 617 

if 
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Figure 16. Experimental set-up utilizing solar concentrator hellostat. 
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Figure 18. Experimental set-up utilizing arc lamp solar 
simulator. 
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B. CESIUM VAPOR CELLS 

The first design selected for the cesium absorption cell is shown in 

Figure 20. This cell was fabricated by ILC Corporation (Sunnyvale, 

California). It has two 1-in. -diameter sapphire windows (brazed to the 

cell with an alkali-resistant braze and an internal depleted uranium 

getter. ILC vacuum baked the cell, loaded it with cesium, and sealed it. 

They specified the cell to be capable of operating to at least 650°C. 

The main problem with this cell was window fogging, which continued 

to occur with this cell design even when high-power window heaters were 

used. Also, the entrance window on the cell cracked twice when the solar 

simulator was used with it, as shown in Figure 21, causing a long delay 

* 

during repair of the cell at ILC Corporation. This cell design ulti- 
mately had to be abandoned because of these window failures, but a 
series of experiments were carried out with the call in the solar con- 
centrator; these are described in Section 3.C. 

We now feel that window fogging in this cell could probably be 
reduced substantially by using an appendage in the cell that is kept 
colder than the cell body, thus serving as a reservoir for liquid cesium. 
Also, since the entrance window broke only when the solar simulator was 
used, we suspect the cause to be the absorption in the sapphire of far-IR 
radiation from the arc lamp. This could be eliminated by using an IR 
blocking filter. 

After the second window failure, we decided to convert the absorp- 
tion cell design to that of a heat pipe oven. Heat pipe ovens offer 
several advantages over equilibrium ovens, such as uniform temperature 

vapor zones, vapor purification, and isolation of the vapor from the 
31 

windows. The latter feature is especially important for our applica- 
tion because it allows using ordinary glass windows and should eliminate 
the window fogging problem by separating the cesium vapor zone from the 

zone adjacent to the windows. Since the vapor pressure in a heat pipe 

31 

oven is determined by the pressure of the inert gas in the oven, the 
cell design included a valve to permit changing the vapor pressure and 
moving the cell outside when the solar concentrator is used. 
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NOTES: 

1. DEVICE VACUUM BAKED AT 600°C 

2. FILLED WITH 10 mg CESIUM 

Figure 20. Original ILC Corp. cesium cell design. 
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Figure 21. Cracked entrance window of ILC cell. 
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The design of Che heac pipe oven underwent three major changes during 
the experiments. This design evolution, shown in Figure 22(a, b, and c) , 
was an effort to eliminate window fogging in the low-aspect-ratio cell 
geometry required to interface converging radiation with the cell. Win- 
dow fogging was finally reduced to a tolerable level by using the conical 
baffles illustrated in Figure 22(c). These baffles allow the incident 
radiation to pass into the central region of the cell relatively unob- 
structed while reducing the size of the convection cells in the end 
region of the cell to a small enough scale so that the cesium is deposited 
on the baffles rather than on the window. Operating the cell in a ver- 
tical configuration (with the lower end of the central zone sealed off) 
reduced the window fogging further and allowed the cesium deposited on 
the baffles to drain back into the central region. 

The cesium vapor pressure could be maintained at temperatures 
between 500 and 600°C for indefinite periods. The window does not fog 
as long as the illumination time by the solar simulator is kept to a 
period of 'vl to 2 sec at a time. The focal volume region in the vapor 
is viewed with an optical multichannel analyzer through the entrance 
window using a carefully aligned, shielded optical train. A photograph 
of the cell in the experimental setup in which the final demonstration 
of plasma heating was achieved is shown in Figure 23. 

The final piece of experimental apparatus fabricated on this pro- 
gram was a heat pipe cell with transverse viewing ports and precision- 
made conical baffles. This cell is shown in Figures 24 and 25. Funding 
limitations prevented carrying out any experiments with this cell. 

C. EXPERIMENTAL RESULTS 

The first experimental indications of plasma formation by solar 
heating were obtained with the equilibrium cesium cell used with the 
Fresnel lens solar concentrator, as shown in Figure 16. The cell was 
operated at 'U500°C, and the cesium dimer absorption bands near 4500 A 

O 

and 5200 A were monitored as the intensity of the incident focused sun- 
light was brought from an attenuated level to full illuminance. The 
results are shown in Figure 26. The reduced absorption gradients as a 
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Figure 23. Experimental set-up which demonstrated plasma heating. 

Solar simulator is at the top, heat pipe cell is ?t 
center, and monocrometer is at the left of the 
photograph. 


49 


Ml 


1 


*743-10 





50 





crc.c J.U PAGE IS 

OF pcca QUALITY 


13300 0740-11 



Figure 25. New heat pipe cell showing conical 
baffles. 
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Figure 26. High-temperature cesium absorption spectra under partial 
illumination from the solar concentrator (a,b) and under 
full ^377-W illumination by the concentrator (c,d). 
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function of wavelength evident at full illuminance indicate thermal 
dissociation of the dimer ground state caused by heating of the incident 
radiation field. Window fogging problems prohibited any quantitative 
interpretation of the data from these experiments. 

When the window fogging problem was finally eliminated with the 

baffled heat pipe cell apparatus shown in Figure 23, three different 

spectral signatures of plasma heating were, observed. An independent 

result for the plasma temperature was obtained from two of these. The 

first spectral signature of plasma heating is the reduction of the 

ground-state dimer absorption caused by thermal dissociation. This is 

shown in Figure 27. In Figure 27(a), a characteristic dimer absorption 

band is seen near 6300 A when the incident solar simulator radiation is 

-3 

attenuated by about a factor of 10 before entering the plasma cell. 
Under illumination by unattenuated radiation from the simulator, the 
Figure 27(b) results show that the absorption strength on these bands 
has been reduced by at least 90%. The associated plasma temperatures can 
be estimated from the absorption strength dependence on dimer concentra- 
tion and the temperature dependence of the equilibrium fraction of 
dimers : 


Y~ = exp (-CJL [Cs 2 ]) 
0 

(26) 

[Cs 2 ] - K [Cs] 2 

(27) 

[Cs] = [ Cs ] Q (T q /T) , 

(28) 


where a is the absorption cross section, L is the optical path length in 
the vapor, T^ is the oven temperature, T is the plasma temperature, K is 
the equilibrium constant, [Cs 2 ] is the ground-state dimer concentration, 

and [Cs] is the equilibrium monomer concentration at T = T . In the 

0 18 -3 0 

experiment, L - 10 cm, T q = 500°C, [Cs] Q = 1 * 10 cm , 

K = 1.3 • 10" 23 exp(5237/T) cm 3 , and 0 (A « 6300 A) = 1 • 10" 16 cm 2 (see 
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Figure 27. 


Plasma formation signatures: dimer dissociation. 
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Figure 1). From the data in Figure 25, I/I Q (6300 A) S 0.9 
within the experimental error. This implies 

T > 1528°K . (29) 

The second spectral signature of plasma heating is shown in 
Figure 28, together with a reference cesium lamp emission spectrum. The 

o o 

absorption lines present in the plasma cell spectrum at 6723 A and 6973 A 
2 2 

are. due to the 7d D + 6p P transitions on the monomer, illustrated in 
Figure 29. Thus, the absorption dips are a direct indication that the 
resonance level of cesium is populated. The absorption cross section for 
these lines is given by 

°*A(^) a ' (30) 

6 

The A coefficient is estimated as A =12*10 , and AX from the spectra is 
M.0 A. These give 


a ~ 2.2 • 10~ 15 cm 2 . (31) 

2 

The population of the 6p P resonance level can be estimated from 

l/l o = exp(-cr[Cs6p 2 P]L) . (32) 

With L = 10 cm, as before, and I/I Q = 0.86, we estimate 

[Cs6p 2 P] a 6.8 • 10 12 cm 3 . (33) 

The temperature, assuming a Boltzmann population distribution (which is 
valid at these densities) , can be estimated from 

[Cs6p 2 P] = 3[Cs o ](T o /T) exp{~1.4/T(eV)l , (34) 
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Figure 28. Plasma formation signatures: excited state absorption. 
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where the factor of 3 is the appropriate degeneracy factor. This gives 
a temperature estimate of 


T s 1320°K (35) 

The third spectral signature of plasma heating is the emission seen _ 

at the 8, 523- A resonance line as well as the emission seen at 8761 A from 
2 2 

the 6d D 6p P transition shown in Figure 30. The latter emission indi- 
cates population of a level at 2.8 eV (see Figure 29). The incomplete 
knowledge of the effective solid angle of the plasma emission restricts 
us to simply stating that the above calculated values of the plasma tem- 
perature are consistent with the emission strengths observed on these 
lines . 

We were able to compare these temperature estimates with the plasma 

temperature predicted by the theoretical model under our experimental 

conditions. Measurements of the incident radiation flux from the solar 

simulator in the focal volume of the baffled heat pipe cell in these 

experiments are shown in Figure 31. The power density is seen to be 

2 

reduced by the baffles to ^20 W/cm , with the focal region degraded to 
about a 2 cm diameter spot. At this power level, dimensional scale and 
oven temperature Cv»500°C) theory predicts that the energy loss from the 
plasma is dominated by radiation. The radiative energy balance results 
of Figure 7 predict a plasma temperature of , ul800 o K under these 
conditions. 
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Figure 30. Plasma formation signatures: line emission. 
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SECTION 4 


CONCLUSIONS 

The theoretical value of 1800°K for the teraperqtura predicted to 
result from heating cesium vapor by the solar simulator represents good 
agreement with the two experimental measurements of the temperature con- 
sidering the difference of the xenon arc lamp spectra from the assumed 
black body form and uncertainties of the effective optical paths through 
the focal volume. This agreement confirms the predictive power of the 
model and establishes increased confidence in our prediction that a 

plasma temperature of n»2800°K will result when the cell is Illuminated 

2 

with sunlight concentrated by the Fresnel lens to 'X'SOO W/cm . 

Thus, solar-driven MHD remains an extremely attractive approach to 

the generation of electric power in space, potentially representing an 

otderS“0 ^magnitude improvement in cost and weight factors compared to 

photovoltaic for large-scale power production. 

There are three well-defined research milestones that now remain 

in order to bring the high-temperature solar MHD converter concept to 

the stage of an operating solar electric generator. The first is an 

experimental confirmation of our theoretical prediction that an 

, u2800 o K plasma temperature can be reached through heating cesium vapor 

2 

by sunlight concentrated to ^300 W/cm by a Fresnel lens. The second 
milestone is the establishment of a theoretical model of the complete 
solar-heated-plasma MHD cycle with which efficiency and weight optimi- 
zation studies can be carried out. Finally, the third milestone is the 
design and fabrication of a prototype high-temperature solar MHD electric 
generator. 
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